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Abstract 

><\ 

Y , Three- and four-jet production is measured in deep-inelastic ep scattering at low x and Q"^ 

Oh| with the HI detector using an integrated luminosity of 44.2 pb^^. Several phase space 

regions are selected for the three-jet analysis in order to study the underlying parton dy- 
namics from global topologies to the more restrictive regions of forward jets close to the 
proton direction. The measurements of cross sections for events with at least three jets are 
^ ' compared to fixed order QCD predictions of 0{af) and 0{a^) and with Monte Carlo sim- 

^ . ulation programs where higher order effects are approximated by parton showers. A good 

\^ I overall description is provided by the 0{a^) calculation. Too few events are predicted at 

CN ' the lowest x ~ 10~^, especially for topologies with two forward jets. This hints to large 

contributions at low x from initial state radiation of gluons close to the proton direction and 
unordered in transverse momentum. The Monte Carlo program in which gluon radiation is 
generated by the colour dipole model gives a good description of both the three- and the 
four-jet data in absolute normalisation and shape. 
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1 Introduction 

The HERA electron-proton collider has extended significantly the available kinematic range 
for tests of Quantum Chromodynamics. The high centre of mass energy of 3 19 GeV allows for 
deep-inelastic scattering (DIS) at a large negative four momentum transfer squared Q^ > 4 GeV^ 
on partons which carry a very small fraction x of the proton momentum down to values of 10"''. 
This is the domain of high parton densities in the proton dominated by gluons and sea quarks. 
In addition, DIS at low x corresponds to scattering at high 'j*p centre of mass energies and is 
therefore intimately linked to the high energy behaviour of QCD. 

Many of the available calculations for DIS processes make use of coUinear factorisation [1]: 
the cross sections are expressed as a convolution of hard partonic subprocesses with proton par- 
ton density functions (PDFs). The latter describe the probabilities to find partons in the proton 
which carry a fraction x of the proton momentum. The separation of the calculation into two 
pieces is specified by the factorisation scale /i^: initial state radiations from the proton with vir- 
tualities above this scale are treated in the hard partonic part while those below are absorbed in 
the PDFs. For inclusive DIS, Q^ provides the natural scale for /ij, i.e. fij = Q"^. The evolution 
of the PDFs with /ij is generally described by the DGLAP [1] equations. To leading logarith- 
mic accuracy this is equivalent to the exchange of a parton cascade, with the exchanged partons 
strongly ordered in virtuality up to Q^. For low x this becomes approximately an ordering in 
kx, the transverse momentum of the partons in the cascade as shown in figure \T\ However, at 
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Figure 1: Exemplary parton cascade diagram for DIS at low x: in the approximation of the 
DGLAP leading log Q^ resummation the emitted gluons are strongly ordered in their transverse 
momenta kx- 



low X the coUinear factorisation scheme may break down. The DGLAP leading logarithmic 
approximation neglects topologies of gluon radiation with unordered kq-. These appear in the 
full perturbative expansion as (log 1/x)" terms which are naturally expected to become large at 



low X. At very low values of x it is believed that the theoretically most appropriate description 
is given by the BFKL evolution equations [2], which resum large logarithms of 1/x. The BFKL 
resummation imposes no restriction on the ordering of the transverse momenta within the parton 
cascade. Compared to the DGLAP approximation more gluons with sizable transverse momen- 
tum are emitted near the proton direction, referred to in the following as the forward direction. 
This should lead to a significantly increased rate of forward jets [3]. A promising approach to 
parton evolution at low and larger values of x is given by the CCFM [4] evolution equation, 
which, because it uses angular-ordered parton emission, is equivalent to the BFKL ansatz for 
X — *> 0, while reproducing the DGLAP equations at large x. 

Higher order calculations in the coUinear factorisation scheme can also improve the treat- 
ment for the low x region, since the log 1/x terms can be treated up to the given order of ag. It is 
an interesting question how well such fixed order approximations can work or if one would still 
need a full log 1/x resummation to all orders as provided by the BFKL approach. For inclusive 
DIS ep — > eX, QCD analyses (see e.g. [5-9]) were performed using coUinear factorisation cal- 
culations in next-to-leading order (NLO) C(as)LJ and/or next-to-next-to-leading (NNLO) O^al) 
for the hard subprocess with parton densities matched to that order. These calculations are able 
to describe the inclusive DIS data from HERA and fixed target experiments over a large range 
in Q^ down to Q^ ~ 2 GeV^ from the largest to the lowest covered x values. 

Final states with jets in DIS are an ideal tool to investigate low x dynamics: the jets can be 
used to tag higher order processes and furthermore provide direct access to the outgoing hard 
partons. The HI and ZEUS measurements of dijet production at low x [10-12] and of inclusive 
forward jet production [13-15] show that the leading order (LO) 0{as) calculations based on 
DGLAP greatly underestimate the data. NLO 0{al) calculations can account for some of the 
LO deficiencies, but the description remains unsatisfactory at low x and Q^. In the present 
paper events with at least three- or four-jets in the final state are investigated. In contrast to 
inclusive jets and dijets, three- and four-jet final states require the radiation of at least one and 
two hard gluons respectively in addition to the qq pair from the dominating hard boson-gluon- 
fusion scattering process 7* (7 -^ qq (see figure |2l). Therefore three- and four-jet processes are 
ideally suited to study the gluon emissions and the underlying parton dynamics in the proton. 
Three-jet cross sections in DIS have been measured previously, both by the HI [16] and the 
ZEUS [17, 18] collaborations. In these analyses the leading jets were required to have a large 
transverse momentum of at least about 7 GeV. All measured cross sections were found to be 
well described by NLO 0{al) predictions in the coUinear factorisation scheme. 

This paper presents a new measurement of three-jet production. The analysis is performed in 
an extended phase space, covering jets with low transverse momenta down to 4 GeV, and based 
on a three times larger luminosity than used in the previous HI publication [16]. The analysis 
reaches values of x as low as x = 10~^. In addition, cross sections for events with at least four 
jets are measured for the first time in DIS. The data are compared with the NLOJET-i-i- [19] 
fixed order calculations in the coUinear factorisation scheme. This program provides LO 0{al) 
predictions for the three-jet case. In addition NLOJET-i-i- is the only available program which 
provides perturbative calculations for jet cross sections in hadronic collisions to 0{al) accu- 
racy. This corresponds to NLO and LO precision for the three- and four-jets cases, respectively. 



'The notation used throughout this paper is that for any calculation labelled 0(a"), the prediction of a cross 
section includes all orders up to n. 



The two Monte Carlo generators RAPGAP [20] and DJANGOH [21], which were able to de- 
scribe reasonably well inclusive forward jet and dijet production at low x, are also tested. 

The sensitivity to deviations from the DGLAP approach may be increased by selecting 
kinematic regions where gluon radiation is suppressed for this approximation. This is the case 
for events with a hard forward jet and a large separation in rapidity to a central parton system. 
Two different subsets of the inclusive three-jet sample are studied: one sample with one forward 
jet and two central jets and one with two forward jets and one central jet. 



2 Kinematics and Measurement Observables 

Figure [2l shows two examples of DIS processes dominating the production of three or more jets 
at low X. The diagrams contribute to order a"^ and a^ to the cross section, respectively. The 
radiated gluons are predominantly emitted in the forward direction whereas the quarks from the 
hard scattering process are mostly central. 





Figure 2: Examples of leading order (left) and next to leading order (right) diagrams for three-jet 
production in DIS at HERA with one and two radiated gluons, respectively. 



The kinematic variables which describe the hard electron-quark scattering process are the 



negative four momentum transfer squared Q^ 



-q = {k — k') of the exchanged virtual 



photon (7*), the Bj0rken variable x = Q'^/{2pq), and the inelasticity y = {qp)/{kp), where 
k, k', p and q denote the four momenta of the incoming and outgoing positron, the incoming 
proton and the exchanged photon, respectively. The three variables are related by Q^ = xys, 
where s denotes the fixed ep centre of mass energy squared. Jets are defined in this analysis in 
the 7*p centre of mass system. The observables used to characterise the jets are their transverse 
momentum p^ in the 7*p centre of mass frame and their pseudorapidity r] in the laboratory 



system. The topology of a three-jet system is fully specified by the following four canonical 
variables [23]: the scaled energy of the jets X; = 2E'J{E[ + E'^ + E'^) (i = 1,2; E[> E'^> E'^) 
and the two three-jet angles 0' and ip' as defined in figure [3l These variables are measured in 
the three-jet centre of mass frame. Four-jet events have additional degrees of freedom. The two 
jets with the lowest dijet mass are combined in order to use the same variables as in the three-jet 
case. 



Three- 
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Figure 3: Definition of the angles 0' and i})' in the three-jet rest system [23]. The 3-vector 
V beam '^^ defined by p '^^^^ = P [n,i - P L,2 where p („ ^ and p (^ ^ are the 3-momenta of 
the two incoming interacting particles in the three-jet centre of mass frame. They are sorted 
with respect to their energy in the laboratory frame: i^in, i > -Ein, 2- The incoming interacting 
particles are the exchanged virtual boson and the parton from the proton side (predominantly 
a gluon). The latter is assumed to move parallel with the proton and to carry a fraction of its 
momentum reconstructed as a;giuon = a; ■ (1 + (s/Q^)), where s denotes the squared centre of 
mass energy of the three-jet system. 



3 QCD Predictions 



The RAPGAP [20] and DJANGOH [21] Monte Carlo event generator programs are used in this 
analysis to estimate the corrections that must be applied to the data for the finite acceptance, 
efficiency and resolution of the detector. The two programs are also used to provide predictions 
that can be compared with the data. Both programs generate hard QCD 2^2 subprocesses 
(e.g. 7* (7 -^ qq) which are convoluted with the CTEQ5L [24] set of parton distributions for the 
proton. The factorisation and renormalisation scales are set to jif = fi^ = Q"^ for DJANGOH 
and /if = /i^ = Q^ + Pr for RAPGAP, where pT is the transverse momentum of the outgoing 
hard partons. RAPGAP includes resolved photon processes using the SaS 2D photon parton 
distribution functions [25], which were found to give a good description of the effective photon 
structure function as measured by HI [26]. Higher order QCD effects that produce further 
hard outgoing partons are generated in both RAPGAP and DJANGOH by parton showers: in 
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RAPGAP the showers are ordered in the transverse momenta (fc-r) of the emissions, according 
to the DGLAP leading log Q^ approximation. DJANGOH uses the colour dipole model (CDM) 
[22], in which partons are generated by colour dipoles, spanned between the partons in the 
cascade. Since the dipoles radiate independently, there is no kx ordering. For the hadronisation, 
the Lund string fragmentation [27] is used both for RAPGAP and DJANGOH. QED radiative 
corrections are applied in DJANGOH using the HERACLES [28] program and are neglected 
for RAPGAP. The DJANGOH predictions are referred to as CDM in the following. 

Fixed order QCD predictions at parton level are calculated using the NLOJET++ [19] pro- 
gram, which is able to predict three-jet cross sections in LO 0{o?^ or NLO 0{c?^ and four-jet 
cross sections in LO 0{p?^. The renormalisation and factorisation scales are set to 



m ^-^ 

i=l 



m ' ^ 



with m = 3 for the three-jet and m = 4 for the four-jet cross sections and A^jet being the number 
of jets fulfilling the applied jet cuts. The value of a^Mz) is fixed to 0. 1 1 8 and the CTEQ6M [7] 
proton parton density parameterisations are used. The NLOJET-I-+ parton level cross sections 
are corrected bin-by-bin for hadronisation effects using the CDM simulation as discussed in 
detail in section |43I Two uncertainties are considered for the NLOJET-i-i- cross sections: The 
uncertainty due to missing higher orders is estimated by recalculating the cross sections with 
the scales fif and jij. varied by a common factor of 2 or 0.5. Hadronisation uncertainties are 
estimated by determining the corrections to the hadron level alternatively with RAPGAP and 
taking 50% of the difference between the corrections from CDM and RAPGAP as systematic 
error. 



4 Experimental Procedure 
4.1 The HI Detector 

A detailed description of the H 1 detector can be found in [29] . Here, a brief account of the com- 
ponents most relevant to the present analysis is given. The HI coordinate system convention 
defines the outgoing proton beam direction as the positive z axis, also referred to as the 'for- 
ward' direction. The polar angle 9 is defined with respect to this direction. The pseudorapidity 
is given by r/ = — In tan(^/2) . 

The central ep interaction region is surrounded by two large concentric drift chambers 
(CJCs), operated inside a 1.16 T solenoidal magnetic field. Charged particles are measured 
in the pseudorapidity range —1.5 < r] < 1.5 with a transverse momentum resolution of 
(^{Pt)/Pt - 0.005 Py / GeV © 0.015. Two additional drift chambers (CIZ, COZ) complement 
the CJCs by precisely measuring the z coordinates of track segments and hence improve the de- 
termination of the polar angle. Multi-wire proportional chambers (MWPC) provide fast signals 
for triggering purposes. 



A finely segmented electromagnetic and hadronic liquid argon calorimeter (LAr) covers the 
range —1.5 < rj < 3.4. The energy resolution is a{E)/E = O.ll/y^E/GeV for electromag- 
netic showers and a{E)/E = O.bO/y^E/GeV for hadrons, as measured in test beams [30]. A 
lead/scintillating fibre calorimeter (SpaCal, [31]) covers the backward region —A<r]< —1.4. 

The data sample of this analysis was collected using a trigger which requires the scat- 
tered positron to be measured in the SpaCal, at least one high transverse momentum track 
(Pt > 800 MeV) to be reconstructed in the central tracking chambers and an event vertex to 
be identified by the MWPCs. The trigger efficiency is higher than 85% for the whole analysis 
phase space. 

The ep luminosity is measured via the Bethe-Heitler Bremsstrahlung process ep -^ ep 7, 
the final state photon being detected in a crystal calorimeter at z = —103 m. 



4.2 Event Selection and Kinematic Reconstruction 

A detailed account of this analysis can be found in [32]. The data used in this analysis were 
taken in the 1999 and 2000 running periods, in which HERA collided protons with an energy of 
920 GeV with 27.5 GeV positrons, corresponding to a centre of mass energy of ^/s = 319 GeV. 
The integrated luminosity of the data is 44.2 pb^^ DIS events are preselected requiring a 
scattered positron measured in the SpaCal with an energy E^. > 9 GeV. The polar angle O^, of 
the scattered positron is determined from the cluster position in the SpaCal and the z position of 
the event vertex reconstructed with the central tracking chambers. The observables y, Q^ and x 
are derived from the electron kinematics 



sin — , Q = 4Ee, beam -Ee COS — , X 



^e, beam -^ -^ 1/S 

where Ee, beam is the positron beam energy. The kinematic range is chosen to be 5 GeV^ < 
g2 < 80 GeV^ 0.1 < y < 0.7, 10"^ < a; < 10'^ ^nd 156° < ^e < 175°. 

The hadronic system, containing the jets, is measured with the LAr and SpaCal calorimeters 
and the central tracking system. Calorimeter cluster energies and track momenta are combined 
using algorithms which avoid double counting [33]. Jets are formed from the hadronic final 
state particles boosted to the 7*p rest frame. The inclusive k^ cluster algorithm [34] is applied 
with a separation parameter of 1.0. The px weighted recombination scheme is used in which 
the jets are treated massless. The jets are ordered with respect to their transverse momentum in 
the 7*p rest frame {p^^ > Pxi+i)- Only jets with a transverse momentum p^- of at least 4 GeV 
and a pseudorapidity in the range —1 < r]i < 2.5 are considered for the analysis. The latter 
cut ensures the jets to lie well within the acceptance of the LAr calorimeter. At least three jets 
are required which fulfil these cuts. It is demanded in addition that p^^ + p'^^ > 9 GeV. The 
applied cuts ensure a good correlation between jets at detector level and hadron or parton level 
and allow for comparison of the data to the NLO 0{al) calculation. In addition to the above 
selections one of the three leading p^ jets has to lie in the central region — 1 < t^j < 1.3. This 
ensures a good trigger efficiency. After all cuts, 38400 events are selected with at least three 
jets. 5900 of these events contain at least four jets. 



Phase Space Definition 



0.1 <y <0.7 

5 GeV^ < g2 < 80 GeV^ 

10-^ <x< 10^2 

156° < ^e < 175° 

K > 9 GeV 



iVjet > 3 



p^. > 4 GeV 

Pri +Pt2>^ GeV 
-1 < ?7j < 2.5 



at least one jet in the range 

-1 <7]i< 1.3 



Table 1: The kinematic domain in which the cross sections are measured. 

4.3 Cross Section Measurement 

The kinematic region for which the cross sections are measured is given in table[TJ All cross sec- 
tions are given as bin-averaged differential cross-sections defined at the level of stable hadrons. 
Therefore the data are corrected for all detector effects, using Monte Carlo simulations. For 
each generated event the response of the HI detector is simulated in detail including trigger 
effects. The events are then subjected to the same reconstruction and analysis programs as the 
data. For each measurement bin a correction factor is calculated as the ratio of simulation en- 
tries at stable hadron level to that at detector level. The same inclusive kr algorithm is applied 
at the hadron and detector levels. The detector correction factors are determined using the CDM 
simulation which is found to give a better description of the jet topologies than RAPGAP. The 
Monte Carlo events are weighted in a few variables to adjust their kinematic distributions to the 
data. These variables are the p^ of the leading jet, r]i — r/2, % + r]2 and Q^. After weighting the 
simulations provide a reasonable description of the shapes of all data distributions. The detector 
correction factors have been studied in detail for all distributions. They vary between 0.6 and 
1.2 for events with at least three-jets (0.4 and 1.2 for events with at least four-jets) and show 
a smooth behaviour. Further small corrections are applied to the data to take QED radiative 
effects into account. The data are corrected to the QED Born level using the CDM simulation. 
A correction factor is determined for each measurement bin separately. 

For comparison with the data, the fixed order NLOJET-i-i- parton level calculations are cor- 
rected to the stable hadron level by application of hadronisation correction factors Chad- These 
corrections are estimated bin-by-bin using the weighted CDM simulation. Jets are obtained at 
the parton level using the inclusive kx algorithm, both in NLOJET-i-i- and CDM. For CDM the 
algorithm is applied to the partons after the parton showering step. As just mentioned, the de- 
tector and hadronisation corrections are calculated using the weighted Monte Carlo simulation 
events. However, the unweighted Monte Carlo predictions are compared to the data, as will be 
shown in section [51 

The correlations between the jets at the different levels have been studied in detail us- 
ing Monte Carlo simulated events. According to CDM, for the phase space given in table [H 
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73-85% of the selected detector level jets can be associated with a hadron level jet within a cone 
A_R = y^A?7^ + A0^) < 0.4 around the detector level jet and with a hadron level transverse 
momentum of p^ > 1.5 GeV. This fraction of "matched" detector jets decreases to 65-75% 
at the parton level. The lowest matching fractions are observed for the more forward jet pseu- 
dorapidities. These migrations dilute the interpretability of the data in terms of the underlying 
partonic dynamics and must be well controlled. The fraction of unmatched jets observed in 
RAPGAP agrees with CDM to better than 30% everywhere. Taking the differences from RAP- 
GAP and CDM as the absolute uncertainty of the number of unmatched jets and assuming this 
number to directly propagate into the measured cross sections a maximal possible cross section 
error of 25% is derived. This possible error will be only considered in the discussion of the 
results (section [5]), whenever a large excess of data over NLO prediction is observed. It is not 
included in the standard determination of systematic cross section errors which is described in 
the following. 



4.3.1 Estimate of Systematic Errors 

The errors of the measured differential cross sections are separated into statistical errors of the 
data (5stat5 additional uncorrelated errors (5uncorr5 accounting for the statistical errors of the Monte 
Carlo samples used to determine the various correction factors, and systematic errors. The latter 
are separated into two contributions: a global normalisation error 5norm and a correlated error 
5corr which affccts only the shape of the cross section distributions. 

The effects of systematic uncertainties on the cross sections are evaluated by applying appro- 
priate variations to the Monte Carlo simulations. The following sources of error are considered: 

• LAr hadronic energy scale: The absolute hadronic energy scale of the LAr calorimeter is 
known to 4% accuracy. This is the dominating uncertainty for the determination of the 
energy of the jets studied in this analysis. 

• SpaCal electromagnetic energy scale: The energy of the scattered positron is known 
within a 2% uncertainty. 

• Positron angle measurement: The uncertainty in the measurement of the polar angle of 
the scattered positron is 1 mrad. 

• Track contribution to combined objects: The uncertainty of this contribution is estimated 
by varying the momenta of all contributing tracks by ±3%. 

• Trigger efficiency: The simulated trigger efficiencies are compared with the efficiencies 
determined from data, using monitor trigger samples. Agreement is found within 3%. 

• Luminosity measurement: The measurement of the integrated luminosity is accurate 
within 1.5%. 

• Radiative correction: The uncertainty of the radiative correction factors is estimated to be 
2% [35]. 
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Error source 


Systematic variation 


"^corr 


'Jnorm 


LAr hadronic energy scale 


±4% 


±3% 


+140/ 
-12/0 


SpaCal em energy scale 


±2% 


±2% 


+0.10/ 
-0.8/0 


Track contribution to comb. obj. 


±3% 


±1% 


±2.5% 


Positron angle 


±1 mrad 


±0.5% 


±1.3% 


Trigger efficiency 


±3% 


< 0.1% 


±3% 


Luminosity 


±1.5% 


< 0.1% 


±1.5% 


Radiative correction 


±2% 


< 0.1% 


±2% 


Model uncertainty 


±50%(CDM-RAPGAP) 


±8% 


n% 


Sum 






tf,% 



Table 2: Sources of correlated systematic uncertainties and the resulting errors on the cross sec- 
tions for events with at least three jets. The first column contains the error source and the second 
the range in which the quantity was varied to account for its uncertainty. The two remaining 
columns give the typical correlated (i^corr) and global normalisation (^norm) uncertainties on the 
cross sections which arise from each error source. 

• Model uncertainty: The cross sections are corrected to hadron level using the weighted 
CDM simulation. The uncertainty of these corrections is estimated by calculating the cor- 
rection factors with the weighted RAPGAP simulation and taking 50% of the difference 
to CDM as systematic uncertainty. 

Typical values for the correlated uncertainties and the global normalisation error on the cross 
sections for events with at least three jets are given in tabled The systematic errors are domi- 
nated by the LAr hadronic energy scale. The second largest contribution stems from the model 
uncertainty. The total global normalisation error is ^14%. The systematic uncertainties for the 
cross sections of events with at least four jets are found to be of similar or somewhat larger size 
as those for events with at least three-jets; the total global normalisation error is ^^9%. 



5 Results 



The measured cross sections are shown in the figures |4] to [121 and listed in the tables [3] tcflOl 



5.1 Cross Sections for Events with at Least Three Jets 



Differential cross sections are presented in figures |4] to [6] for events with three or more jets 
as a function of the number of jets (A'jet), the Bj0rken variable x, the pseudorapidities of the 
three jets and the variables characterising the topology of the three jets in the three-jet centre 
of mass frame (scaled jet energies and three-jet angles). The kinematic range for which the 
cross sections are determined is specified in table [H The figures also show the predictions 
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of the NLOJET++ fixed order QCD prediction in LO C(a2) and NLO 0{a^), corrected for 
hadronisation effects. The theory error, including scale variations and hadronisation correction 
uncertainties added in quadrature, exceed the measurement uncertainty. Figure |4] (left) shows 
the distribution of the number of jets found in the selected events which extends up to A^jct = 6. 
For this distribution the predictions of the NLO 0{a^) calculation and of the two Monte Carlo 
programs RAPGAP and CDM are shown. The NLO 0{ag) prediction agrees for A'jet = 3 and 
underestimates the rate of events with 4 jets by a factor 2.6. It does not produce any events with 
more than four jets as expected. 

The differential cross sections shown in figures |4] to [6] are not described by the LO C(ag) 
QCD predictions neither in shape nor in magnitude. The main discrepancies are seen at low 
X and for forward jets (large positive t]) where far too few events are predicted. The NLO 
0{al) prediction improves the situation considerably in all regions where deviations from LO 
are observed. A similar improvement was already noticed in the previous HI three-jet analysis 
[16], which was restricted to the phase space of large invariant three-jet masses above 25 GeV. 
For that phase space the 0{al) calculation could still describe the three-jet data down to the 
smallest x = 10^"^. However, for the present analysis without such mass cuts applied, at very 
small X < 2 ■ 10^"^ the calculation undershoots the data, which lies approximately at the upper 
edge of the total theoretical error band, by a factor of about 0.6. In the complementary region 
X > 2 ■ 10~^, the description is reasonable; this was also observed in the recent ZEUS multi- 
jet analysis [18] which was restricted to this phase space. In summary, a large deficit of the 
NLO prediction persists only at low x and for forward jets. This is the kinematic region where 
unordered gluon radiation is expected to enhance the jet production [3]. The shapes of the 
cross section distributions for the three-jet topological variables shown in figure [6] are all well 
described by the NLO 0{a^) prediction, only the global normalisation of the calculation is 
somewhat too low. 



5.2 Forward Jet Subsample 

The observed excess of data versus QCD predictions in the region of forward jet rapidities and 
at low X is investigated here in further detail. The sample of selected events with at least three 
jets fulfilling the criteria presented in table[T]is reduced to a subsample by requiring that at least 
one of the three leading jets be forward and carry a large proton momentum fraction: 

E- t 
r]> 1.73 and Sjet = ^ > 0.035 . 

-^p, beam 

Here Ej^t is the jet energy in the laboratory frame and Ep^ beam = 920 GeV the proton beam 
energy. Further requirements are applied to obtain two separate subsamples. In the sample with 
one forward jet and two central jets (f+2c) the other two leading jets are required to lie in the 
central pseudorapidity range —1 < ?7jct < 1- In the sample with two forward jets and one central 
jet {2f+c) it is demanded that one of the three leading jets is a central jet with — 1 < r^jet < 1 and 
the remaining leading jet must be in the more forward region r^jct > 1 (for this second forward 
jet no cut is applied on a^jot)- The fraction of jets due to gluon radiation is expected to be larger 
for forward jets than for central jets. This is confirmed by a study of the parton composition 
of three-jets in the CDM simulation. Therefore the f+2c sample will have many events with 

13 



only a single radiated gluon (as for the left diagram in figure O while the 2f+c selection will 
have a larger fraction of events with two radiated gluons (as for the right diagram in figure [2l). 
Cross section measurements as a function of x, % and p^^ for the /+2c and 2f+c samples are 
presented in figures |7] and [8l The additional global normalisation errors of tu:% for the /+2c 
sample and li5% for the 2f+c selection are not shown in the figures. The fixed order NLO 
0{a^) prediction gives a rather good description of the /+2c sample. The step from LO 0{a'^) 
to NLO 0{a^) improves the agreement mainly at very low x < 2 ■ 10^'', where a remaining 
deficiency of ~30% is observed. Towards larger x values both the LO and the NLO calculation 
fall off somewhat less steeply than the data and are too high for the largest covered x values 
X > 2 ■ 10^^. For the 2f+c selection an even more dramatic change is observed at low x from 
C(ag) to 0{a^): the discrepancy at x < 2 ■ 10""^ is reduced from a factor of 10 to 2.6. The 
large remaining deficiency exceeds the combined error of prediction and data and is thus highly 
significant. It can also not be explained by a possible additional maximal cross section error of 
25% due to detector jets which cannot be matched with hard partons (as discussed in the last 
paragraph of section |43I) . This data excess provides a strong hint for missing higher order QCD 
corrections, i.e. beyond 0{a^), in this forward gluon radiation dominated phase space. Note 
that for processes with two radiated gluons, the 0{al) calculation can only provide a leading 
order perturbative estimate. 

Excesses which are probably related to the one reported here were observed in the forward 
jet analyses from HI [13] and ZEUS [15]. In these analyses the topologies of three jets were 
investigated for events containing a dijet system in addition to a forward jet. The 0{a^) pre- 
dictions were found to undershoot the data in the region where all three jets tend to go forward. 
However, in these analyses the data were either integrated over a larger x range or restricted 
to somewhat larger x values, which might explain why the excesses are less prominent than 
observed in the present measurement. 



5.3 JMonte Carlo Program Predictions for Events with at Least Three Jets 

The cross sections for events with at least three jets are compared to predictions from RAPGAP 
and CDM. The jet multiplicity shown in figure |4] (left) is described well by CDM while RAP- 
GAP falls off too steeply. The overall normalisation of the RAPGAP and CDM predictions is 
found to be too low by 55% and 5% respectively. In the following the Monte Carlo predictions 
are normalised to the total measured cross section in order to compare only the shapes of the 
cross sections. Figures |9] and \T0\ show the comparison as a function of the Bj0rken scaling 
variable x, the difference of the pseudorapidity of the two leading p^ jets (Ar/ = r/i — r]2), the 
variables p^^, X2 and the two three-jet angles cos^' and cos^^'. RAPGAP fails to describe the 
X and Ai] distributions. CDM on the other hand gives a very good description of almost all 
observables besides p^-^^, where it predicts too many high momentum jets (p^^ > 15 GeV). The 
three-jet angular distribution cos ip' is also described rather poorly by both RAPGAP and CDM. 
A separate check of the cross sections for p^^ > 20 GeV reveals significant deviations to CDM 
for the shapes of various distributions, especially for x and r]i, as presented in figure [TT] In this 
domain RAPGAP describes these distributions well. 



14 



5.4 Cross Sections for Events with at Least Four Jets 

A subsample of events with four or more jets is also studied. All selection criteria of the three- 
jet sample listed in table[T]have to be fulfilled. In addition, at least one more jet has to be found 
which satisfies the standard jet cuts p^ > 4 GeV and — 1 < r/ < 2.5. As already mentioned the 
NLOJET++ 0{a^) calculation can only provide a LO prediction for the final state with four jets 
which is by far too low as can be seen in figure |4] (left). Thus in the following the comparisons 
of the measured four jet cross sections are restricted to the CDM and RAPGAP predictions, 
where parton showers approximate higher orders and can lead to large jet multiplicities. The 
total cross section predicted by CDM for events with four or more jets agrees well with the data 
while RAPGAP is too low by a factor of ~ 2.9, as shown in figure |4] (left). Differential cross 
sections as a function of p^^, rfi — 774, X2 and cos 6' for events with at least four jets are shown 
in figure [T2I and compared to the predictions by the two Monte Carlo generators normalised to 
the data. RAPGAP fails to describe the shapes of the differential distributions, again with the 
exception of the momentum distributions of the jets. CDM on the other hand disagrees with the 
data in the p^ distributions but describes the scaled energies of the four jets correctly. It also 
gives a very good description of all other distributions. 



6 Summary 

This paper presents a new measurement of three-jet production in DIS at low x and Q^. The 
measurement is carried out in an extended kinematic phase space covering lower jet transverse 
momenta compared to previous three-jet analyses. Very small x values are reached down to 
X = 10~^. The first measurement of four-jet production in DIS is also presented. Three- and 
four-jet final states require the radiation of at least one respectively two hard gluons from the 
initial state proton, in addition to the qq pair from the dominating hard boson-gluon-fusion 
scattering process 'j*g — > qq, and are therefore well suited to study parton dynamics at small x. 

The measurements are compared with the NLOJET-i-i- [19] fixed order QCD calculations. 
A remarkable result of the present analysis is the success of the next-to-leading order 0{a^) 
calculation for the cross sections of events with at least three jets. The inclusion of diagrams 
with two radiated gluons improves dramatically the agreement with the data compared to the 
0{a^) prediction which is far too low especially at small x. A similar improvement was already 
noticed in the previous HI three-jet analysis [16]. In the present analysis, extending to lower 
invariant three-jet masses, an excess is observed of the data compared to the 0{a^) prediction 
at the lowest x ~ 10^^. This excess is found to be enhanced and to become highly significant 
for topologies with two forward jets and one central jet. Excesses which are probably related, 
albeit less significant, were observed in the forward jet analyses from HI [13] and ZEUS [15]. 
The new analysis corroborates the hypothesis that the DGLAP leading log Q^ approximation 
starts to break down in the region of the lowest accessible x ~ 10'^ for Q^ > A GeV^, at least 
up to the order 0{al) for which calculations are presently available. In other words in a sizable 
fraction of events, which is much larger than predicted, two or more gluons are radiated from 
the initial state proton which are unordered in their transverse momentum, i.e. they all have 
relatively large transverse momenta. For events with at least four jets the 0{al) prediction is 
also too low, as expected, since the calculation can provide only a leading order estimate. 
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The new data presented here are also compared with two Monte Carlo simulation programs 
implementing hard QCD 2— >2 processes complemented by parton showers modelling higher 
order effects. The RAPGAP [20] program, using kx ordered parton showers and including 
resolved photon processes, fails to describe the data. On the other hand the DJANGO [21] pro- 
gram with non fc^-ordered gluon radiation as implemented in the colour dipole model (CDM) 
gives a remarkably good description of the measured cross sections for events with at least three 
and four jets and even for events with higher jet multiplicities. The remaining discrepancies at 
high p^ require further studies. 

The three- and four-jet production is investigated further by a detailed study of the jet topolo- 
gies as represented by three-jet angles and scaled momenta. The best description for the case of 
three-jet production is obtained by the 0{al) NLOJET-i-i- calculation, significantly better than 
that provided by CDM. 
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Figure 4: Differential cross sections as a function of the number of jets A'jet found in the 
events and the Bj0rken scaling variable x. The results are obtained for the selected events 
with at least three jets in the kinematic range listed in table [U The cross sections are bin- 
averaged and plotted at the respective bin centers. The inner error bars represent the statistical 
error of the data, the total error bars correspond to the statistical and uncorrected systematic 
uncertainties added in quadrature. The hatched error bands show the estimate of the correlated 
systematic uncertainties. The data have an additional overall global normalisation error of '^11% 
(not shown). The dark shaded (inner) error band shows the NLO 0{al) prediction with the 
uncertainty due to the hadronisation corrections, the light shaded (outer) band includes the 
scale uncertainty added in quadrature. The dashed line represents the LO C(«s) prediction. 
The latter is not shown in the A'jct distribution which instead is compared with the two Monte 
Carlo programs RAPGAP (direct + resolved, dashed line, labelled as 'RG d+r') and CDM (solid 
line). 
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Figure 5: Differential cross sections as a function of the pseudorapidity r]i for each of the three 
leading jets {p^^ > p*rp^ > p^g). The data are compared to the LO 0{al) and to the NLO 0{al) 
predictions. See the caption of figure |4] for further details. 
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Figure 6: Differential cross sections as a function of the scaled energies X[ and Xg of the two 
leading jets (E[ > E'2 > E!^ in the three-jet centre of mass frame) and the two angles 9' and ^' 
as defined in figure [3l The data are compared to the LO 0{a'^) and the NLO 0{a^) predictions. 
See the caption of figure |4]for further details. 
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Figure 7: Differential cross sections as a function of the Bj0rken scaling variable x for the two 
selected subsamples of events with at least three jets: the sample with one forward jet and two 
central jets (f+2c) on the left, the sample with two forward jets and one central jet (2f+c) on 
the right. The additional global normalisation errors of the data (1}4% for the /+2c sample and 
'til% for the 2f+c selection) are not displayed. See the caption of figure|4]for further details. 
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Figure 8: Differential cross sections as a function of the pseudorapidity r]i and the transverse 
momentum p^^ of the leading jet for the two selected subsamples of events with at least three 
jets: the sample with one forward jet and two central jets (f+2c) on the left, the sample with two 
forward jets and one central jet (2f+c) on the right. The additional global normalisation errors 
of the data (^^4% for the /+2c sample and ^^5% for the 2f+c selection) are not displayed. See 
the caption of figure |4] for further details. 
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Figure 9: Differential cross sections as a function of the Bj0rken scaling variable x and the 
difference of the pseudorapidities of the two leading jets for the selected events with at least 
three jets in the kinematic range listed in table \\\ The cross sections are bin-averaged and 
plotted at the respective bin centers. The inner error bars represent the statistical error of the 
data, the total error bars correspond to the statistical and uncorrelated systematic errors added 
in quadrature. The correlated systematic errors are shown by the hatched error band. The 



additional global normalisation error (; 



-16c 
14' 



of the data is not displayed. The data are compared 



to predictions from DJANGO (CDM) (solid line) and RAPGAP (dashed line). Both Monte 
Carlo predictions are scaled to match the total data cross section: CDM by a factor of 1.05 and 
RAPGAP by 1.55. 
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Figure 10: Differential cross sections for the selected events with at least three jets as a function 
of the transverse momentum of the leading jet p^^, the scaled energy X2 of the next to leading 
jet in the three-jet centre of mass frame and the two angles 9' and tp' as defined in figure |3l The 
data are compared to the predictions from the Monte Carlo programs RAPGAP and CDM. Both 
Monte Carlo predictions are scaled to match the total data cross section: CDM by a factor of 
1.05 and RAPGAP by 1.55. See the caption of figure|9]for further details. 
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Figure 1 1 : Differential cross sections as a function of the Bj0rken scaling variable x and the 
leading jet pseudorapidity rji for the selected subsample of the events with at least three jets: 
the leading jet is required to have a transverse momentum p^^ above 20 GeV. The data are 
compared to the RAPGAP and CDM predictions. Both Monte Carlo predictions are scaled to 
match the total data cross section: CDM by a factor of 0.41 and RAPGAP by 1.95. See the 
caption of figure |9] for further details. 
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Figure 12: Differential cross sections for events with at least four jets as a function of the 
transverse momentum p^^ of the leading jet, the pseudorapidity difference of the leading and 
the fourth jet r]i — 774 (with p^-^ > p^^ > p^g > p^^), the scaled energy X2 and the angle 9' 
as defined in figure [3l For the determination of X2 and 9' the four jets are reduced to three by 
combining the two jets with the lowest dijet mass. The resulting three jets (labelled with ") are 
then sorted with respect to their energy in the three-jet centre of mass frame {E" > E'^ > E'^) 
and both variables are calculated as for the three-jet case, for instance X2 = p„ J, p„ . For the 



details of the plotting of the data see the caption of figure[9j The additional global normalisation 
error of the data (^^19%) is not displayed. The data are compared to the RAPGAP and CDM 
predictions. Both Monte Carlo predictions are normalised to the total cross section of the data. 
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A Cross Section Measurement Tables 



Cross sections for events with at least three jets 



Njet 


da/dNjet 


Ostat 


'^uncoiT 


OcoiT 


Chad 




(pb) 


(%) 


(%) 


(%) 




3 


4.8- 10^ 


1 


1 


+2 
-3 


0.71 ±0.07 


4 


5.6-101 


2 


3 


+ 15 
-14 


0.65 ±0.03 


5 


6.0 


4 


10 


+9 
-6 


0.70 ±0.10 


6 


0.25 


11 


36 


+55 
-52 


0.64 ±0.30 


X 


da/dx 


Ostat 


f^uncorr 


Ocorr 


Chad 




(pb) 


(%) 


(%) 


(%) 




[0.0001,0.0002) 


1.3-106 


2 


3 


+ 10 
-12 


0.74 ±0.08 


[0.0002,0.0005) 


7.1 - 10^ 


1 


2 


+3 
-4 


0.70 ±0.08 


[0.0005,0.001) 


2.3- 10^ 


1 


2 


+5 
-4 


0.69 ±0.07 


[0.001,0.002) 


6.2 - 10'' 


2 


2 


+ 12 
-10 


0.66 ±0.07 


[0.002,0.005] 


7.7-103 


2 


3 


+23 
-16 


0.65 ±0.04 



Table 3: Bin-averaged differential cross sections for the selected events with at least three jets 
in the kinematic range listed in table [T] The cross sections are defined at the hadron level and 
are given as a function of the jet multiplicity A^jet and the Bj0rken scaling variable x. The 
following cross section uncertainties are indicated: statistical (^stat)^ uncorrected systematic 



-16/ 



and correlated systematic (5corr)- The additional global normalisation uncertainty of 
_il% is not included in the table. The correction factors Chad for the effect of hadronisation 
and the associated uncertainty are also given. They are applied to the NLOJET++ parton level 
calculations. 
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Cross sections for events with at least three jets 



m 


da/diji 
(pb) 


Ostat 
(%) 


^uncorr 

(%) 


Ocorr 

(%) 


Chad 


[1.0,-0.5) 
[-0.5,0.0) 
[0.0,0.5) 
[0.5,1.0) 
[1.0,1.5) 
[1.5,2.0) 
[2.0,2.5] 


78 
150 
191 
204 
187 
136 
139 


2 
2 
1 
1 
2 
2 
2 


3 
3 
2 
2 
2 
3 
3 


+11 

-8 
+9 
-7 
+2 
-2 

+3 
-3 

+2 
-3 

+4 
-5 
+4 
-5 


0.64 ±0.09 
0.69 ±0.08 
0.70 ±0.09 
0.74 ±0.06 
0.70 ±0.07 
0.69 ±0.05 
0.72 ±0.07 


m 


dcr/d?72 
(pb) 


Ostat 
(%) 


f^uncorr 

(%) 


Ocorr 
(%) 


Chad 


[1.0,-0.5) 
[-0.5,0.0) 
[0.0,0.5) 
[0.5,1.0) 
[1.0,1.5) 
[1.5,2.0) 
[2.0,2.5] 


105 

157 
188 
193 
179 
128 
130 


2 
2 
1 
1 
1 
2 
2 


3 
3 
2 
2 
3 
3 
3 


+4 
-4 
+6 
-7 
+4 
-5 
+3 
-4 
+4 
-3 
+6 
-4 
+8 
-6 


0.62 ±0.08 
0.67 ±0.09 
0.73 ±0.09 
0.72 ±0.05 
0.72 ±0.06 
0.70 ±0.07 
0.72 ±0.07 


m 


da/drjg, 
(pb) 


Ostat 
(%) 


^uncorr 

(%) 


Ocorr 
(%) 


Chad 


[1.0,-0.5) 
[-0.5,0.0) 
[0.0,0.5) 
[0.5,1.0) 
[1.0,1.5) 
[1.5,2.0) 
[2.0,2.5] 


81 

122 

177 

186 

188 

159 

174 


2 
2 
2 
1 
1 
1 
1 


4 
3 
3 
3 
2 
3 
2 


+5 
-4 
+4 
-3 
+3 
-4 
+ 1 
-1 
+4 
-3 
+2 
-1 
+3 
-4 


0.60 ±0.09 
0.67 ±0.06 
0.71 ±0.06 
0.74 ±0.08 
0.74 ±0.19 
0.69 ±0.08 
0.72 ±0.04 



Table 4: Bin-averaged differential cross sections for the selected events with at least three jets 
in the kinematic range listed in table[I] The cross sections are defined at the hadron level and are 
given as a function of the pseudorapidities 7]^^ 4=1,2,3 of the three leading jets in the lab frame. The 
following cross section uncertainties are indicated: statistical ((5stat), uncorrelated systematic 
^uncorr and Correlated systematic (^corr)- The additional global normalisation uncertainty of 
^^11% is not included in the table. The correction factors Chad for the effect of hadronisation 
and the associated uncertainty are also given. They are applied to the NLOJET++ parton level 
calculations. 
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Cross sections for events with at least three jets 



X[ 


da/d{X[) 


Ostat 


^uncorr 


Ocorr 


Chad 




(pb) 


(%) 


(%) 


(%) 




[0.6,0.7) 


667 


2 


3 


+3 
-1 


0.59 ±0.08 


[0.7,0.8) 


1940 


1 


2 


+ 1 
-3 


0.65 ±0.08 


[0.8,0.9) 


2650 


1 


2 


+ 1 
-1 


0.72 ±0.06 


[0.9,1.0] 


1260 


1 


2 


+5 
-4 


0.83 ±0.06 


^2 


da/d(X^) 


Ostat 


^uncorr 


^corr 


Chad 




(pb) 


(%) 


(%) 


(%) 




[0.5,0.6) 


726 


2 


3 


+2 
-1 


0.78 ±0.08 


[0.6,0.7) 


2570 


2 


2 


+3 
-3 


0.68 ±0.07 


[0.7,0.8) 


1710 


1 


2 


+3 
-3 


0.69 ±0.07 


[0.8,0.9) 


397 


2 


4 


+9 
-6 


0.74 ±0.04 


[0.9,1.0] 


11 


12 


15 


+ 13 
-12 


0.70 ±0.06 


cos 6*' 


dcr/d(cos6'') 


Ostat 


f^uncorr 


^corr 


Chad 




(pb) 


(%) 


(%) 


(%) 




[-1.0,-0.8) 


359 


2 


3 


+5 
-3 


0.70 ±0.06 


[-0.8,-0.4) 


307 


2 


2 


+7 
-6 


0.71 ±0.09 


[-0.4,0.0) 


216 


1 


3 


+6 

-7 


0.71 ±0.07 


[0.0,0.4) 


220 


1 


3 


+ 10 
-11 


0.74 ±0.04 


[0.4,0.8) 


288 


1 


3 


+2 
-3 


0.72 ±0.08 


[0.8,1.0] 


300 


2 


3 


+6 
-5 


0.60 ±0.08 


cosip' 


dcr/d(cosV'') 


Ostat 


^uncorr 


^corr 


Chad 




(pb) 


(%) 


(%) 


(%) 




[-1.0,-0.8) 


437 


2 


2 


+6 
-7 


0.70 ±0.06 


[-0.8,-0.4) 


240 


2 


3 


+4 
-3 


0.71 ±0.06 


[-0.4,0.0) 


221 


1 


3 


+3 
-3 


0.72 ±0.08 


[0.0,0.4) 


220 


1 


3 


+3 
_2 


0.69 ±0.07 


[0.4,0.8) 


236 


2 


3 


+6 
-5 


0.70 ±0.08 


[0.8,1.0] 


441 


1 


2 


+7 
-8 


0.70 ±0.07 



Table 5: Bin-averaged differential cross sections for the selected events with at least three jets 
in the kinematic range listed in table[T] The cross sections are defined at the hadron level and are 
given as a function of the three-jet topological observables X[, X'^, cos 9' and cos ip' as defined 
in section [2l The following cross section uncertainties are indicated: statistical (4tat), uncor- 
rected systematic (5uncorr and correlated systematic (^corr)- The additional global normalisation 
uncertainty of li4% is not included in the table. The correction factors Chad for the effect of 
hadronisation and the associated uncertainty are also given. They are applied to the NLOJET-i-i- 
parton level calculations. 
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Cross sections for events with at least three jets 



Pti 


dc^/d(PTi) 


Ostat 


^uncoiT 


dcorr 


Chad 


(GeV) 


(GcV^^ • pb) 


(%) 


(%) 


(%) 




[5,8) 


60 


1 


2 


+2 
-2 


0.67 ±0.09 


[8,11) 


56 


1 


2 


+2 

— 2 


0.69 ±0.07 


[11,15) 


29 


1 


2 


+2 
— 2 


0.73 ±0.06 


[15,20) 


9.5 


2 


3 


+8 
-8 


0.77 ±0.05 


[20,25) 


2.6 


4 


5 


+8 
-10 


0.79 ±0.01 


[25,30) 


0.82 


7 


7 


+ 15 
-10 


0.80 ±0.02 


[30,45) 


0.18 


9 


8 


+ 13 
-13 


0.81 ±0.01 


[45, 60] 


0.016 


32 


17 


+39 
-26 


0.80 ±0.09 


m - m 


dcr/d(?7i - 772) 


Ostat 


^uncorr 


f^corr 


Chad 




(pb) 


(%) 


(%) 


(%) 




[-3.0,-2.3) 


25 


3 


6 


+27 
-16 


0.70 ±0.05 


[-2.3,-1.5) 


66 


2 


4 


+2(1 

-14 


0.72 ±0.07 


[-1.5,-0.8) 


113 


1 


3 


+9 
-8 


0.72 ±0.09 


[-0.8,0.0) 


144 


1 


2 


+5 
-7 


0.71 ±0.07 


[-0.0,0.8) 


146 


1 


2 


+ 11 
-12 


0.69 ±0.07 


[0.8,1.5) 


118 


1 


3 


+5 
-8 


0.70 ±0.05 


[1.5,2.3) 


76 


3 


3 


+2 
-3 


0.68 ±0.08 


[2.3,3.0] 


30 


3 


5 


+ 11 
-9 


0.70 ±0.05 



Table 6: Bin-averaged differential cross sections for the selected events with at least three jets 
in the kinematic range listed in table[T] The cross sections are defined at the hadron level and are 
given as a function of the transverse momentum p^^ of the leading jet in the 7*j9 centre of mass 
frame and the pseudorapidity difference r/i — 7/2 of the two leading jets in the lab frame. The 
following cross section uncertainties are indicated: statistical (^stat)^ uncorrected systematic 
^uncorr and Correlated systematic (^corr)- The additional global normalisation uncertainty of 
t)[i% is not included in the table. The correction factors Chad for the effect of hadronisation 
and the associated uncertainty are also given. They are applied to the NLOJET++ parton level 
calculations. 
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Cross sections for events with one forward jet and two central jets 



X 


dcr/da; 
(pb) 


"stat 
(%) 


^uncorr 

(%) 


•^corr 

(%) 


Chad 


[0.0001,0.0002) 

[0.0002,0.0005) 

[0.0005,0.001) 

[0.001,0.002) 

[0.002,0.005] 


1.0 -lO^ 

6.0 • lO'^ 

2.1 • 10"* 
4.9 • 10^ 
4.0 • 10^ 


5 
3 

5 
7 
12 


9 
6 

7 
9 
16 


+7 
-5 
+2 
-5 
+3 
-2 
+ 17 
-13 
+ 14 
-17 


0.82 ±0.10 

0.72 ±0.07 
0.69 ±0.04 
0.66 ±0.08 
0.74 ±0.10 


m 


dcr/d?7i 
(pb) 


Ostat 
(%) 


^uncorr 

(%) 


Ocorr 
(%) 


Chad 


[-0.5,0.0) 
[0.0,0.5) 
[0.5,1.0) 
[1.0,1.5) 
[1.5,2.0) 
[2.0,2.5] 


8.4 
12.9 
16.9 
0.0 
10.4 
39.0 


7 
6 
5 

6 
3 


13 
10 
8 


10 
6 


+ 14 
-10 
+24 
-15 
+5 
-8 
+0 
-0 
+6 
-7 
+4 
-7 


0.76 ±0.07 
0.77 ±0.20 
0.76 ±0.03 
0.00 ±0.00 
0.78 ±0.04 
0.72 ±0.08 


Pti 
(GeV) 


dc^/d(PTi) 
(GeV~^ • pb) 


Ostat 
(%) 


fJuiicorr 

(%) 


Ocorr 
(%) 


Chad 


[5,8) 
[8,11) 
[11,15) 
[15,20) 
[20,25) 
[25,30) 
[30, 45] 


1.9 
4.4 
3.9 
1.4 
0.36 
0.11 
0.027 


6 
4 
4 
6 
11 
21 
29 


12 
7 
7 
9 

13 
19 
22 


+ 10 
-8 
+4 
-6 
+ 14 
-2 
+ 14 
-17 
+29 
-34 
+ 18 
-15 
+28 
-16 


0.70 ±0.20 
0.71 ±0.05 
0.74 ±0.07 
0.76 ±0.09 
0.82 ±0.01 
0.87 ±0.10 
0.82 ±0.07 



Table 7: Bin-averaged differential cross sections for the selected events with at least three jets 
in the kinematic range listed in table [TJ In addition one of the three leading jets is required 
to be a forward jet with ^jct < 20° and to carry a large fraction of the proton beam energy 
a; Jet > 0.035. The two other jets are required to lie in the central region of the HI detector 
—1 < r^jct < 1 (sample with one forward jet and two central jets). The cross sections are 
defined at the hadron level and are given as a function of the Bj0rken scaling variable x and the 
leading jet observables: pseudorapidity rfi in the lab frame and transverse momentum p^^ in 
the 7*p centre of mass frame. The following cross section uncertainties are indicated: statistical 
((^stat), uncorrelated systematic (5uncorr and correlated systematic {Scott)- The additional global 
normalisation uncertainty of ^{4% is not included in the table. The correction factors Chad for 
the effect of hadronisation and the associated uncertainty are also given. They are applied to the 
NLOJET++ parton level calculations. 
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Cross sections for events with two forward jets and one central jet 



X 


dcr/dx 
(pb) 


"stat 
(%) 


*^uncorr 

(%) 


'^corr 

(%) 


Chad 


[0.0001,0.0002) 

[0.0002,0.0005) 

[0.0005,0.001) 

[0.001,0.002) 

[0.002,0.005] 


1.3 
8.9 
3.0 
9.3 
1.6 


10^ 
lO'^ 
10-^ 
103 
103 


4 
3 
3 
4 
6 


8 
5 
5 
6 

8 


+9 
-7 
+7 
-9 
+8 
-6 

+ 11 
-10 

+ 13 
-13 


0.76 ±0.09 
0.77 ±0.03 
0.72 ±0.05 
0.70 ±0.06 
0.69 ±0.04 


m 


(pb) 


Ostat 
(%) 


^uncorr 

(%) 


^corr 
(%) 


Chad 


[1.0,-0.5) 
[-0.5,0.0) 
[0.0,0.5) 
[0.5,1.0) 
[1.0,1.5) 
[1.5,2.0) 
[2.0,2.5] 


2.9 
7.7 
9.3 
10.4 
15.3 
30.0 
62.6 


11 

7 
6 
5 
5 
4 
2 


20 
14 
12 
11 
8 
6 
4 


+29 
-22 
+6 
-9 

+27 
-17 
+ 10 
-9 
+6 
-5 
+5 
-7 
+3 
-4 


0.70 ±0.01 
0.82 ±0.01 
0.92 ±0.01 
0.88 ±0.04 
0.72 ±0.10 
0.71 ±0.04 
0.73 ±0.05 


Pti 
(GeV) 


dc^/d(PTi) 
(GeV"^ • pb) 


Ostat 
(%) 


fJuiicorr 

(%) 


Ocorr 
(%) 


Chad 


[5,8) 
[8,11) 
[11,15) 
[15,20) 
[20,25) 
[25,30) 
[30, 45] 
[45, 60] 


2.8 
6.4 
5.4 
2.4 
0.8 
0.28 
0.065 
0.006 


4 
3 
3 

4 
7 
11 
14 
60 


10 
6 
5 
6 
9 
11 
11 
22 


+ 10 
-10 

+4 
-4 

+3 
-5 
+8 
_2 

+8 
-7 

+5 

-11 

+ 11 
-10 

+57 
-32 


0.72 ±0.05 
0.68 ±0.05 
0.76 ±0.06 
0.78 ±0.04 
0.82 ±0.09 
0.79 ±0.01 
0.78 ±0.01 
0.86 ±0.30 



Table 8: Bin-averaged differential cross sections for the selected events with at least three jets 
in the kinematic range listed in table [IJ In addition one of the three leading jets is required 
to be a forward jet with 6'jct < 20° and to carry a large fraction of the proton beam energy 
xjct > 0.035. One of the two other jets has to fulfil r^jct > 1 and the other —1 < ?7jet < 1 
(sample with two forward jets and one central jet). The cross sections are defined at the hadron 
level and are given as a function of the Bj0rken scaling variable x and the leading jet observ- 
ables: pseudorapidity rji in the lab frame and transverse momentum p^^ in the 7*p centre of 
mass frame. The following cross section uncertainties are indicated: statistical (5stat), uncor- 
rected systematic ^uncorr and correlated systematic ((5corr)- The additional global normalisation 
uncertainty of ^14% is not included in the table. The correction factors Chad for the effect of 
hadronisation and the associated uncertainty are also given. They are applied to the NLOJET-i-i- 
parton level calculations. 
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Cross sections for events with at least three jets and p^^ > 20 GeV 



m 


da/di]i 
(pb) 


Ostat 
(%) 


^uncorr 

(%) 


(%) 


Chad 


[-0.5,0.0) 
[0.0,0.5) 
[0.5,1.0) 
[1.0,1.5) 
[1.5,2.0) 
[2.0,2.5] 


1.4 
4.6 
8.2 
10.3 
9.0 
7.1 


18 
10 

7 
6 

7 
7 


33 
13 

8 
7 
7 
8 


+30 
-15 
+ 13 
-11 
+7 
-7 
+7 
-7 
+ 10 
-8 
+7 
-7 


0.71 ±0.40 

0.88 ±0.40 
0.76 ±0.40 
0.80 ±0.40 
0.80 ±0.40 
0.79 ±0.40 


X 


da/dx 


Ostat 


fJuiicorr 


^corr 


Chad 


[0.0001,0.0002) 

[0.0002,0.0005) 

[0.0005,0.001) 

[0.001,0.002) 

[0.002,0.005] 


(pb) 
4.6 • lO'^ 

2.5 • lO'^ 

9.6 • 103 
2.4- 10^ 
2.0-102 


(%) 

7 
5 
7 
9 
17 


(%) 

7 
5 
6 

10 
20 


(%) 

+8 
-4 
+2 
-3 
+3 
-7 
+8 
-6 
+23 
-26 


0.82 ±0.40 
0.77 ±0.40 
0.81 ±0.40 
0.81 ±0.40 
0.83 ±0.40 



Table 9: Bin-averaged differential cross sections at the hadron level as a function of the Bj0rken 
scaling variable x and the pseudorapidity r]i of the leading jet in the lab frame, The cross sec- 
tions are measured in the kinematic range listed in table[IJ In addition the leading jet is required 
to have a large transverse momentum p^^ > 20 GeV in the 'j*p centre of mass frame. The 
following cross section uncertainties are given: statistical ((5stat), uncorrected systematic ^uncorr 
and correlated systematic (^corr)- The additional global normalisation uncertainty of ]'^}4% is not 
included in the table. The correction factors c^ad for the effect of hadronisation and the associ- 
ated uncertainty are also given. They are applied to the NLOJET-I-+ parton level calculations. 
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Cross sections for events with at least four jets 



VI - m 


dcr/d(77i - ?;2) 


fJstat 


'^uncoiT 


dcorr 




(pb) 


(%) 


(%) 


(%) 


[-3.0,-2.3) 


2.8 


7 


15 


+5 
-12 


[-2.3,-1.5) 


9.5 


4 


9 


+ 10 
-7 


[-1.5,-0.8) 


15.8 


3 


7 


+4 
-6 


[-0.8,0.0) 


14.9 


8 


7 


+2 
-5 


[-0.0,0.8) 


14.2 


3 


7 


+8 
-9 


[0.8,1.5) 


15.0 


3 


7 


+ 1 
-2 


[1.5,2.3) 


7.5 


5 


10 


+ 15 
-10 


[2.3,3.0] 


3.8 


7 


15 


+47 
-19 


X'2 


dCT/d(X^) 


dstat 


'^uncoiT 


Ocorr 




(pb) 


(%) 


(%) 


(%) 


[0.5,0.6) 


127 


3 


7 


+5 
-8 


[0.6,0.7) 


305 


4 


4 


+ 7 
-6 


[0.7,0.8) 


169 


3 


6 


+ 7 
-10 


[0.8,0.9] 


29 


6 


13 


+20 
-12 


cos6»' 


dCT/d(cos6'') 


Ostat 


'^uncoiT 


Ocorr 




(pb) 


(%) 


(%) 


(%) 


[-1.0,-0.8) 


61.7 


10 


8 


+3 
-2 


[-0.8,-0.3) 


25.5 


3 


6 


+ 10 
-13 


[-0.3,0.3) 


23.2 


3 


6 


+ 7 
-5 


[0.3,0.8) 


29.6 


3 


6 


+4 
-3 


[0.8,1.0] 


46.3 


3 


7 


+9 
-9 


Pti 


da/d{p*^^) 


Ostat 


f^uncoir 


Ocorr 


(GeV) 


(GcV^^ • pb) 


(%) 


(%) 


(%) 


[5,8) 


5.6 


5 


7 


+3 
-3 


[8,11) 


6.3 


2 


6 


+4 
-3 


[11,15) 


3.7 


3 


6 


+4 
-6 


[15,20) 


1.4 


4 


8 


+6 
-8 


[20,25) 


0.41 


8 


11 


+ 11 
-8 


[25, 30) 


0.22 


13 


16 


+ 19 
-16 


[30, 45] 


0.028 


20 


17 


+47 
-24 



Table 10: Bin-averaged differential cross sections at the hadron level as a function of the pseu- 
dorapidity difference rji — 774 of the leading and the fourth jet, the topological observables X'2, 
cos 0' and the leading jet transverse momentum p^^ in the 7*p centre of mass frame. The cross 
sections are measured in the kinematic range listed in table [IJ In addition it is required that a 
fourth jet is found in the events which fulfils the standard jet selection cuts p^ > 4 GeV and 
— 1 < r/ < 2.5. The following cross section uncertainties are given: statistical (5stat), uncor- 
rected systematic ^uncorr and correlated systematic (^corr)- The additional global normalisation 
uncertainty of li9% is not included in the table. 
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